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Dyring the present reporting period work was finished on extending 
the Hadley Cell model to investigate the syitinetric baroclinic instability 
of the model. This work was written in the form of a paper which will 
appear in the June 1982 issue of The Journal of The Atmospheric Sciences. 

A copy of the preprint of this paper Is enclosed as Appendix. At the 
present work is being continued to investigate the baroclinic instability 
of an inclined wave for the same Hadley pell model. 

.Also during the reporting period work has begun on the coding of a 
program to solve the nonlinear baroclinic instability of an Eady model. 

The basic expression technique have been established and all the theoretical 
background for this work has been completed. It is anticipated the trail 
runs of the program should begin in the coming months. Initially we will 
investigate the nonlinear evolution of a wave of a single zonal wave number 
and the influence of its hOrmonics. The nonlinear interaction of more 
than one zonal wave will be investigated subsequently. 
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• Abstract 

A stability analysis of a thin horizontal rotating fluid 
layer which is subjected to arbitrary horizontal and vertical 
temperature gradients is presented. The basic state is a non- 
linear Hadley cell which contains both EXtnan and thermal boundary 

layers; it is given in closed form. The stability analysis is 

» 

based on the linearized Navier-Stokes equations, and zonally 
symmetric perturbations in the form of waves propagating in the 
meridional direction are considered. Numerical methods were used 
for the stability .problem. The objective of this investigation 
was to extend previous work on symmetric baroclinic instability 
with a more realistic model. Hence ^ the- study deals with flows 
for which the Richardson number- (based on ’temperature and flow 
gradients at mid-depth) is of order unity and less. The com- 
putations cover ranges of Prandtl number 0.2 a ^ 5, Rossby 

-2 2 -4 -1 

number 10 < z <_10 and Ekman number 10 ^E <10 .It was found, in 

agreement with previous work, that the instability sets in when 
•the Richardson number is close to unity, and that the critical 
Richardson number is a non-monotonic function of the Prandtl 
number. Further, it was found’ that the critical Richardson number 
decreases with increasing Ekman number, until a critical value 
of the Ekman number is reached beyond which the fluid is stable. 

The principal of exchange of stability was not assumed and growth, 
rates were calculated. A wavelength of maximum growth rate was 
found. For our model overstability was not found. Some com- 
putations were performed for Richardson numbers less than zero. 






No discontinuities in growth rates are noticeable when the 
Richardson number changes sign. This result indicates a 
smooth transition from symmetric baroclinic instability to a 
convective instability. 







1. Introduction 

It is well known that statically stable baroclinic flow 
can be destablished by two different mechanisms. The first 
occurs in what is usually known as ordinary baroclinic in- 
stability in which the perturbations of maximum growth rate 
have a wavelength in the zonal direction and the gravest 

i 

possible structure in the meridional direction. (Charney, 

1947; Eady, 1949) . This instability is well understood and 
has been observed experimentally (Hide and Mason, 1975 ). The 
second occurs in what is known as symmetric baroclinic in- 
stability. In this case, the instability mechanism is such 
that the perturbations of maximum growth rate have meridional 
structure but no zonal structure. This study is concerned with 
the latter form of instability. 

Two of the earliest analyses of symmetric baroclinic instability 
are those of Solberg (1936) and Kuo (1956) . However, the first 
definitive work is that of Stone. In a sequence of papers. 

Stone (1966, 1970, 1971) examined the stability of a parallel 
baroclinic flow with respect to inf inite-simal , three-dimensional 
perturbations. Constant vertical shear and contstant temperature 
gradients were assumed and thermal and viscous diffusion effects 
were neglected. Stone found that when the Richardson number, 

Ri, is between 0.95 and 0.25, waves whose axes are in the zonal 
direction have the maximum growth rates. Stone et al . (1969) • * 

and Hadlock et al. (1972) attempted to observe symmetric bf^roclinic 

> n 

instability in the laboratory, and some evidence of the predicted 
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meridional structure was seen but a comprehensive experimental 
study of this instability has yet to be performed. 

McIntyre (1970) performed an analysis of symmetric 
baroclinic instability which included both vertical and 
horizontal shears and allowed for the Prandtl number # d. The 
model considered was infinite in both the horizontal and vertical 
directions. McIntyre found, in the^ absence of horizontal shear, 
that the maximum value of Ri at which symmetric instability sets 
in (the critical Richardson number, Ri ) is a function of cr'and 
can be larger than 0.95. McIntyre also found that the most un- 
stable wavelength is infinite and that both monotonic and 
oscillatory instability can occur. In a related analysis, 

Walton (1975) considered only vertical shear in a finite depth 
fluid. When weak viscous effects are included, Walton found 
that Ri^ is lower for both the monotonic and oscillatory ’ in- 
stabilities. He also found that the most rapidly growing modes 
display a weak functional dependence on the diffusion coefficients 
when those are asymptotically small. 

Recently, Emanuel (1979) presented a, detailed study of 
the linear stability of symmetric perturbations in a baroclinic 
flow with both horizontal and vertical shears. Retaining both 
the viscous and thermal dif fusivities for a vertically bounded 
flow, Emanuel studied: i) hydrostatic disturbances and ii) non- 

hydrostatic disturbances in a neutrally stratified fluid. He 
found that for both cases the critical Richardson number depended 
oh both the thermal and the viscous dif fusivities . He also 

2 
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established the fact that the imost unstable normal mode is 
determined primarily by the depth of the unstable domain and 
the slope of the isentropic surfaces rather than by the 
diffusive properties of the fluid. Weber (1980) in a some- 
what similar analysis# examined the symmetric instability of 
a basic state with a constant vertical shear in the presence 
of horizontal boundaries while taking full account of the 
viscous effects for arbitrary • vertical (stable) stratification. 
Although, he made the simplifying assximption of replacing the 
vertical viscous and diffusive terms in the perturbation 
equations by a wave mode, his results compare favorably with 
those of Emanuel. It should be noted that in both of the above 

studies the analyses were restrictive; Emanuel investigated only 

♦ 

neutral modes by invoking the principal of exchange of stability 
and Weber used an approximate solution to the perturbation 
equations. Emanuel's paper contains an excellent survey of 
previous work on this problem. 

Previous workers have postulated- symmetric baroclinic 
instability as the cause of certain physical phenomena. Stone 
(1967) conjectured that the banded structure of Jupiter's 
atmosphere is due to. this instability. Bennets and Hoskins 
(1979) attributed the origin of rain bands and squall lines to 
this instability. Weber (1980) suggested that the generation of 
roll vortices in the oceans and the atmosphere is due to symmetric 
baroclinic instability. 
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In ail of the theoretical studies cited above, the basic 

state whose stability was investigated was a postulated state 

in which both the ahears and the temperature gradients in the 

vertical and the horizontal directions were taken to be constants 

throughout the fluid. In an experimental apparatus whose upper 

and lower boundaries are stationary, such a basic state is a 

good approximation in the interior of the fluid only; it is not 

a good approximation near the boundaries, especially when the 

Ekman number is not vanishingly small. In the study, presented. 

in this paper, a fluid contained between two horizontal plates 

of infinite extent is considered and the influence of the Ekman 

and the thermal layers on symmetric baroclinic instability is 

investigated. The analysis uses a realistic basic state which is 

obtained through the solution of the governing nonlinear equations of 

motion for abitrary external forcings, e.g., the Ekmi.n number^ 

the vertical stratification, etc. This basic state was used 

in a recent paper by the authors in an investigation of ordinary 
« . 

» 

baroclinic instability (see Antar and Fowlis, 1981) . 

It is obvious that if a direct analog of a laboratory ex- 
periment is required, then this model .is not sufficient since it 
does not allow for the side walls and the side layers. In a 
recent 'study Quon (1980, 1981) simulated laboratory conditions 
for the onset of symmetric baroclinic instability through the 
numerical solution of the full Navier-Stokes , energy equations 
for .the flow in a cylindrical differentially heated, rotating 
annulus. Although such an analysis is of great value in. modeling 




4 
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the true experimental conditions, the amount of computer time 
required for the calculation of each case precludes a thorough 
probing of the parameter space for the criteria for the onset 
of the instability. It is our opinion that the present study 
fills a needed gap between the constant shear, constant tem- 
perature gradients analyses discussed earlier and that of the 
exact solution of Quon. 

The basic state model is reviewed in Section 2 . The 
perturbation equations for the symmetric instability and their 
method of solution are presented in Section 3, In Section 4 
we discuss the results cf the present analysis and 
how they compare with the previous work. Our conclusions 
are presented in Section 5. 


5 




2* The Basic State 

The basic state used in the present stability analysis 
is identical to that used in the work of Antar and Fowlis (1981) 
(hereinafter referred to as I) . This model ia described only 
briefly here; for further details and justifications, the . 
reader is referred to I. We consider a Boussinesq fluid confined 
between two horizontal plates which are set a distance, d, apart. 
The coordinate system used ia rectangular Cartesian with axes 
(x,y,z) corresponding to the eastward, northward and vertical- 
directions, respectively, and with the origin midway between the 
plates. A sketch of the model is shown in Fig. 1. The plates 
are assumed to extend to infinity in the x-direction and to 
large distances in the y-direction. The plates and the fluid 
are taken to rotate as a whole about the vertical axis with a 
constant angular velocity, fi. A temperature gradient in v;hich 
the temperature decreases linearly in the y-direccion is main- 
tained on both plates. To assure that the vertical stratification, 
and hence the Richardson number, is arbitrary and externally * 
fixed, the temperatures of the- upper and lower plates are set to 
. differ uniformly by a constant amount AT for all (x,y), 

* The velocity and temperature fields are governed by the 
• « 

Navier-Stokes, en&r,fy and mass conservation equations. Under the 
conditions assumed here for the basic, two-dimensional, steady 
state, these equations, in a rotating reference frame, reduce to 
the following set; 
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-2V * E3^U/3z^, 

(1) 

2U « ES'^V/az^ - 3p/3y/ 

(2) 

V3T/3y * (E/cc) 8*T/3zS 

(3) 

T * 3p/3z, 

(4) 


where is the velocity vector (U,V)< p is the pressure^ T is 
the !s:emperature and 
e « agy/S^ 

I 

E « v/ndS 

a *= v/k‘, 

are a thermal Rossby number, the Ekman number and the Prandtl 

number, respectively. .In Egs. (l)-{4), length, time, velocity 

and temperature have been made dimortsionless using d, , 

t gd v/fi, and yd, where a is the coefficient of thermal expansion, 

* 

y is the imposed horizontal temperature gradient and g is gravity. 
The main assumptions used in arriving at Egs. (l)-(4) are that 
.the vertical velocity component, W, is negligible and that 

I 

3T/9y *=-Y throughout the region of interest, (see I) . 

The solution of the system (l)-(4) which is consistent 
with the no slip and perfectly conducting conditions at the solid 
boundaries, i.e. , 

U~V«0, T=+ AT/2-y 0 z ® + 1/2, (5) 

* » 

is given by 

•U (z) * - f (z)/8 t z/2 , • (6) 


% 
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V(z) 


g(z)/8, 


(7) 


T(y,z) * -y + zAT + eol2z - f(z))/8, (8) 

where 

f(z) * (ooshR(z+l/2) cosR(z-l/2) - coshR(z-l/2)cosR(z+l/2) ] /h (R) , 

g(z) « Isir,hR(z+l/2) sinR(z-l/2) - sinhR(z-l/2)sinR{z+l/2) ]/h(R) / 

h(R) = 8inh*’R/2) + sin^(R/2), 

R = 

The velocity and temperature profiles given by (6) -(8) are 
shown in Fig. 2 for different values E, e, a, and AT. The 
limits of applicability of this model are discussed .in I. 


8 
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3, The Perturbation Equations 

The instability pro*;^lem is tackled by imposing perturba- 
tions on the basic state presented above. Since we are concerned 
with symmetric perturbations, the perturbation functions will be 
proportional to exp [ii(y - ct) ] which represents a wave with 
meridional wavenumber i propagating in the meridional direction with 
speed c. Introducing the perturbations into the Navier-Stokes, 
energy and mass conservation 'equations and upon linearizing the 
perturbation equations we obtain the following equations for the 
perturbation functions . 

E(D*-i.^)^w + iA[(c-eV) (D^-£,^)W + EwD^Vl -2iJ!,Du-X^Q = 0, (9) 

E(D^ 7 JI^) u + iA (c-eV) u-ewDU 7 ^ Dw « 0, (10) 

E(D*-Jl^) 0 + iAa (c-eV) e- ^CDw - eawDT = 0, (11) 

iiu' + Dw = 0, (12) 

where 
D = d/dz. 

In the above the perturbation velocity is = (u,v,w) and 0 

is -the perturbation temperature. Note that the perturbation 
♦ * 
equations are not hydrostatic . 

To complete the specification of the problem, the following 

boundary conditions are imposed on the perturbation functions. 

u«w=Dw=0=O§z-+ 1/2. (13) 


9 
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The perturbation problem as posed constiitutes an eigen- 
value problem* We take 1 to be real and c » + ic. 

r 1 

to be complex^ where denotes the propagation speed and 
the growth rate. Monotonic instability occurs for ^ 0 while 
Cj, * 0, whereas overstability occurs for c^ * 0 while c^ 0. 

It is tempting at this stage to invoke the principle of exchange 
of stability and to set c = c. =0 and hence look for neutral 
modes alone and reduce the number of parameters involved. Emanuel 
(1979) adopted this procedure. However, since McIntyre (1970) 
had established that for his model (see Section 1) oscillatory 
instability exists, it was decided to leave c^ and as eigen- 
values of the problem. Using this method the analysis will allow 
us to compute growth rates and to determine if overstability 
•exists. 

The solution to the problem (9) - (13) is obtained through 
numerical integration of the full eighth- order differential 
system. This procedure was necessary because the coefficients 
in the equations (i.e., the basic state) are functions of the* 
height, z. If we wish to account fully for the Ekman and thermal 
layers, the solution cannot be obtained in a closed form. Note, 
that if the basic state profiles had been taken to be linear in z, 
as was done in most of the previous work on the symmetric baroclinic 
instability problem r then the set of governing equations would 
reduce to a differential system with constant coefficients whose 
solution could. be determined in a closed form. 
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The numerical solution was obtained using the shooting 
technique. The code was identical to the one used in I# and 
more details of the numerical work are given in I. Briefly, the 
shooting technique was enhanced with the orthonormalization 
method, allowing the stability analysis to .be performed even 
when the Ekman number is very small. A Newton -Raph son method 

t 

was implemented for convergence on the eigenvalues, and an eighth- 
order Runge-Kutta method was used for the integration of the 
initial value problem. 




^ • Results and Discussion 

The results which follow constitute the solution of the 
eigenvalue problem posed by equations (9) ~ (11) with the 
boundary conditions (13) for prescribed values of the 
parameters. The parameters of the problem which appear ex- 
plicitly in the equations are, e, the thermal Rossby number 
E, the Ekman number, a, the Prandtl number, AT, the imposed 
vertical temperature difference, and i and c, the wavenumber 
and wave speed, respectively. Note that the first four 
parameters appear in the basic state as well as the perturbation 
equations, but the last two belong to the stability problem 
alone . 

As was discussed in Section 1, the Richardson number plays 
a significant role in. the symmetric baroclinic instability problem 
While this number does not appear explicitly in our formulation 
of the basic state, it may be defined using the parameters which 
do appear explicitly. The Richardson number is defined by 


Ri 


ag 3 T*/ 32 * 

m*Az*p 


(14) 


where 3T*/3z* and 3U*/3z* are the- dimensional basic state tem- 
perature and zonal flow gradients, respectively. It is obvious 
from this definition and from expressions (6) - (8) that for the 
present model Ri is not a .constant but is a function of z. We • 
may define a "buDc” Richardson number by the value of Ri at 
mid-depth between the horizontal plates, viz. 
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3T*/3z* 

Rifa » ag ( ^ 


Z at 0 • 


In all of the subsequent c’lscussion, Rij^ will be used instead 
of Ri and the subscript b will be dropped* We shall see that 
this definition is not only convenient, but it has dynamical sig- 
nificance in this problem. Substituting the temperature and 
velocity gradients for the basic state# expressions (6) and 
(8) f respectively, into the above definition we obtain, 



AT + ae. 

4 ^ 


(15) 


which is a constant and defined exclusively by the basic state. 

The above definition shows that Ri is a function of the 
explicit parameters of the problem* This relationship .was used 
when choosing values of e, AT and a for each eigenvalue search 
in order to be sure that we were within the correct range of Ri 
for instability* Stone (1970) and McIntyre (1970) showed that 
symmetric baroclinic instability occurs for Ri close to unity 
(see Section 1) . The above definition indicates that for a 
given value of a, there is a range pf. e and of AT for which .the 
critical value of Ri is achieved* 


For 0=1, expression (15) specifies that AT must be less 
than zero to obtain the desired range of values of Ri for in- 
stability regardless of the value of e. However, setting AT< 0 
arbitrarily could lead to an unstable vertical stratification 
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in the interior and hence a negative value of Ri. Instabilities 
occuring under the above conditions with negative Ri are not 
classed as symmetric baroclinic instabilities and are outside 
the range of interest of this paper. Such instabilities 
have been discussed by Hathawayet al . (1979 ,1980). To obtain symmetivic 
baroclinic instability for a«l, AT must be set negative with the 
additional constraint of )4 AT/e| <_!». Clearly, there exists 
a range of values AT which satisfies these two criteria. We decided 
to proceed with AT ~ -0.2 for a * 1, Fig. 3 shows that this 
value of at keeps the interior vertical stratification stable. 

The above arguments reveal that the correct choice of AT for 
Epfflietric instability is a delicate matter and depends strongly 
on a. Expression (15) should be of great value in guiding ex- 
perimental searches for symmetric baroclinic instability in the 
laboratory. 

Fig. 4 shows the growth rates for different values of the 

-3 

meridional wavenumber, i,‘ as a function of e for E = 10 , AT « 

-0.2 and a *= 1. Note the existence of unstable modes over a 
finite range of e. This feature resembles the growth rate 
properties of ordinary baroclinic instability. The unstable 
region in e is bounded by two stable regions, one in which Ri is 
close to unity and the other in which Ri is small or negative 
(Ri *= 1 - 0.8/e) . 

Fig. 4 also shows that as e varies, the growth rate has a 
maximum value for a specific value of the wavenumber. Fig. 5, 
which -is for the same values of the parameters as the results in 
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Fig. 4, is a plot of the maximum growth rate for each wavenumber 
versus wavenumber. The greatest growth rate occurs for about 
* 10. This result is at variance with the findings of Stone 
(1970 ) , who found that the maximum growth rate is achieved as 
i -► and McIntyre (1970) , who found that it is achieved as 
i ■+■ 0. However, since both of these previous studies are for 
essentially inviscid flow, it is concluded that including the 
viscous effects determines the wavelength of maximum growth rate 
for the symmetric baroclinic instability problem. 

Figure 6, which is for the same values of the parameters as 
the results in Pigs. 4 and 5, shows neutral stability curves in 
the (E - e) nvamber plane for three different wavenumbers.’ A sig- 
nificant feature of these curves is the existence of the maximum 
(critical) value’ of E‘ above which' all modes are stable*. Another 
feature is the existence of tiwo stable regions, or two branches of 
the neutral stability curves, for all values of E smaller than the 
mode critical values. Note that these curves are loci of points 
for which c^ = 0. For all the modes shown in this figure -and 
those in the previous and subsequent figures, the value of the wave 
speed, c^, was always found to be zero, indicating that these are 
monotonic instability modes and tha» the principle of exchange of 
stability holds for this problem. .Throughout this work we did 
not find any oscillatory unstable modes. Thus, it is doubtful if 
overstability exists for the present problem. 

The curves in the upper portion of Fig. 6 indicate -that 
E « as E 0 for all of the wavenumbers shown. To observe this 
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criterion in terms of Ri, the neutral curves of Fig. 6 are plotted 
again in the (E-Ri) plane in Fig. 7. Fig. 7 shows that as E 0, 
the upper branches of the curves approach Ri = 1 for all wave- 
lengths and that the region of instability increases with in- 
creasing wavenumber. Weber (1980 ) , in his model for symmetric 
baroclinic instability (see Section 1 ) , found the critical value 
of Ri, Ri«/ to be a monotonically decreasing function of E as the 

w 

latter increases from zero. The results in Fig. 7 show that this 

fundamenti.\l dependence is correct up to the value of E^, but above 

that value there is a region of absolute stability regardless of 

the value of Ri. We believe that the existence of E_ is due to 

c 

the accurate inclusion of the viscous* and thermal diffusive effects 
in our model. None of the previous wor)c included both of 
these effects properly. 

TO investigate the effects of the Prandtl number on the 

stability criteria discussed above, we proceeded to find solutions 

for a « 5. As was discussed earlier, changing the value of a 

necessitates a search for the apprppriate value of the vertical 

temperature difference, AT, such that Ri is in the proper range 

for symmetric baroclinic instability,. , The value chosen for o = 5 

was AT = -4.5, which allows Ri to be in the vicinity of unity 

(Ri - 5 - 18/e) . Fig. 8 'shows the growth rates as a function of 

-3 

£ for the above values of a and AT and for E = 10 . Again, it 

is seen that the unstable modes separate two stables regions 
in ,E (or Ri) . Fig. 8 shows that the maximum growth rate increases 
with increasing wavenumber, Jl. In fact, the total range of 
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instability increases with H. Note that the upper bound of 
instability in e first increases \tfith increasing wave numbers 
up to f « 4 and then decreases, reaching an asymptotic value of 
e * 4,55 for 5. > 6 which corresponds to Ri « 1 for the parameters 
used. Another feature, which is different from the a » 1 case, 
is that the maximum growth rates occur for approximately the same 
value of e for all the wavenumbers shown. For this case this 
value of E translates to a value of Ri = 4i 

Fig. 9, which is for the same values of the parameters as 
the results in Fig. 8, shows the variation of the maximum growth 
rates as a function of the wavenumber, X,. This figure, reveals 
that there is still a wavenumber of greatest growth rate which 
occurs about a * 13.5. However, the decrease in the growth 
rate is much smaller for wavenumbers above the maximum value 

than it is for the a *= 1 results (see Fig. 5) , 

The neutral stability curves in the (E - Ri) plane are 
shown in Fig. 10 for f = 1 and 4. In this figure only the positive 
part of the Ri range is shown. Note that there are two values of 
Ri as E -*■ 0, each corresponding to one wavenumber, and that the 
range of instability is much larger for Sj - 4 than that for f = 1. 

To appreciate more fully the influence of the Prandtl 
number on the symmetric instabiity criteria. Fig. 11 which shows 
the neutral stability curves for £ = 4 in the (E - Ri) plane for 

a *= 1 ( AT = -0.2) , a = 2 (AT = -1.5) and a = 5 (AT = -4.5) was . 
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prepared. Except for a « 1/ only the upper branches of the 

curves appear for Ri > 0. All the curves shown tend to specific 

asymptotic values of Ri as E 0. For o * 1, the asymptotic 

o 

value of Ri is close to 1, while for a 2 and 5, the asymptotic 

w 

values are close to 1.075 and 1.695 respectively. The latter 
two values are significantly greater than i. Thus, the trend 
for the asymptotic value of Ri is to increase with increasing 

C 

* 

0 for a > 1. This fundamental dependence is in agreement with 

work of McIntyre (1970), although his model is somewhat different 

from the model being considered here (See Section 1) . McIntyre 

found that the asymptotic value of Ri is given by 

c 



(1 + 0 ) ^ 
4cr 


(16) 


For a * 1, 2, and 5> this relationship yields values of Ri of 1, 

• . o . 

1.125 and 1 . 8, respectively . These values are close to the values 
found from Pig . 11 . 

Pig. 11 also shows that a:s a increases, the wiiole region 
of instability is markedly increased. The value of E increases 

w 

9 

drastically as 0 is increased but there is no observat.le difference 

in the values of E^ for 0=2 and 5. 

c 

Next, we present results fOi 0 < 1. Again, the values for 
AT were chosen such that the range of Ri includes unity and values 
less than unity. Expression (15) indicates that when c < 1 and 
AT < 0, the maximum possible value for Ri is 0 and thus the range ■ 

does not include unity. However, when 0 < 1 and AT > 0, values of 

• * 

Ri can be obtained which are both greater and less than unity.' We 
proceeded with AT > 0 . 
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Fig. 12 shows the neutral stability curves in the (E - Ri) 

plane for a * 0.5 (AT « 0.2) and o » 0.2 (AT » 0.5). The curves 

are for I » 12 , the value which yielded the maximum value for 

Ri when all the other parameters were held fixed. Again, the 

values for Ri^ are seen to increase to specific asymptotic values 

as E 4* 0. Note that the tendency to an asymptote is not as evident 

for a *= 0.5 as it is for a « O.2.— The main features of Pig. 12 

are the decrease in the value of Ri_ and the decrease in the total 

c 

instability area as o increases to unity. The results also in- 
dicate an increase in with increasing a. 

-4 

Pig. 12 gives the values of Ri^ at E *= 10 for a * 0.5 

and 0.2 as 0.946 and 1.21, respectively. Pig. 13 is a comparison 

-4 

of all our results for the largest asymptotic (E « 10 ) values 

of Ri with McIntyre's results as a function of a. The dots are 
o 

our results and the continuous curve is a plot of expression (16). 

Note that McIntyre's criterion yields a symmetric dependence of 

Ri as a function of a about a * 1. The results of the present 
c 

study agree well with McIntyre for a ^ 1 but give lower values 
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5 . Conclusions 

The stability of a nonlinear rotating Hadley cell has 
been examined with respect to symmetric perturbations. This 
Hadley cell was derived previously {see I) and contains both 
Ekman and thermal boundary layers. The stability analysis is 
based on the Navier-Stokes equations and considers infinitesimal 
waves with vertical and meridional structure but no zonal structure. 
The stability analysis was performed numerically. For examini. q 

■ the parametric behavior of symmetric baroclinic instability/ the 

work is an improvement over all previous work since a realistic 

basic state is used and since the stability analysis includes 

» 

■both viscous and thermal diffusion effects. 

The computations cover the following ranges of the non=- 

-o 2 -4 ■ * -1 

dimensional parameters:. 0.-2 ^ o <5, 10 e <, 10 ,10 1 ^ 10 , 

-4.5 ^ AT 0.5/ and -100 ;< Ri _< 2. It was found that the instability 

sets in when is close to unity and that the critical Richardson 

■ nximber/ Ri / is a strong function of both o and E for the cases 

c 

considered. We extrapolated from these cases to obtain values of 

* ♦ 

Ri\, for the inviscid limit, E 0 / and found that Ri_ has a maximum 
value close to unity for cr * 1.' The general trend of these results 
is in agreement with McIntyre (1970) and for ^ ^ 1 the agreement is 
very good but fore < 1 it is not so good. We should not be suprised 

to find differences between the present and previous work because 
of the differences in the basic states and other aspects of the* 
models’. 
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We also found for fixed 0 that Ri„ decreases with in- 

c 

creasing E until a critical value of E is reached beyond which 

tlie fluid is stable. The principle of exchange of stability was 

not assiuned and growth rates were calculated. A wavelength of 

maximum growth rate occuaing for small wavelengths was found. 

For our model overstability was not found. Some computations were 

performed for Ri < 0 . In particular, for o « 5 the growth rate 

m:iximum occurs at about Ri » -4 . No discontinuities in the growth 
» 

rates are noticeable when Ri changes sigc>. This result indicates 
a smooth transition, from symmetric baroclinic instability to a con- 
vective instability. 


« 
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Legend 

1, A sHetch of the geometry of the model and the 
coordihate system. 

2a. The zonal velocity component of the basic state, 

U, as a function of z, for three values of E. 

2b. The meridional component of the basic state, V, 
as a function of z, for the same values of E as 

, Ji . 

the ^results shown in ' Figure 2a. 

2c. Basic state temperature profiles, T at y *= 0 for 
the same values of E as the results shown in 
Pig . 2a and for o = 1, e. - 1 and AT *= 0. 

3. Basic state temperature profiles, T, at y = 0 for 

three values of AT and for a « 1, e * 1, and 

E * 5 X lO"^. 

4. Growth rates, as a function of c for different 

-3 

values of •£, and for a = 1, AT = -0.2 and E « 10 

The upper abscissa shows the values of Ri corresponding 

to e . 

5. The maximum growth rates as a function of the 
meridional wave number. Z, for the same conditions 
as Fig. 4. 

6. Neutral stability curves, c^ « 0, as functions of 

E and e for different values of i and for a * 1 and 
AT *= -0.2. 

7. The neutral stability curves of Fig. 6 as functions . 
of E and Ri. 
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Figure 8. 


Figure 9. 
Figure 10. 


Figure 12 . 


Figure 13. 


Growth rates, as a function of e for different 

values of I and for a * 5, AT •= -4 . 5 and E * 10 , 

The upper abscissa shows the values of Ri corresponding 
to e. . 

The maximum growth rates as a function of I for the 
same conditions as Fig. 8. 

Neutral stability curves as functions of E and Ri 
for different values of Z and for a = 5 and AT * -4.5. 
Neutral stability curves as functions of E and Ri for 
different values of o 1) and for JL = 4. The 
dashed lines are McIntyre's results. 

Neutral stability curves as functions of E and Ri for 

different values of a (< 1) and for iZ, = 12. The 

dashed line is McIntyre's result. 

A comparison between the critical value of the 

Richardson number, Ri_» as calculated from Eq. (16) 
j c 

( ) and the results presented in this paper 

-4 

as a function of a for E =.10 (0). 


25 




ORIGINAL PAGE BS 
OF POOR QUALITY 


k 


% 1 

V * 

y/y//////////////// 


; 

: 

1 

y — 

r 

U‘ \ 

y // Tj -TV - AT /2 y ////////////////////a 
L 


FIGURE 1 


a 

i 

s 


n 

u. 

H 


ii 


26 


f ORIGINAL PAGE m 

OF POOR QUALnYj 



27 










ORIGfrJAl PAGS tS 
OF POOR QUALiTir 



FIGURE 2c 


29 


S 


































University of Tennessee 
Space Institute 
Tullahoma, Tennessee 37388 

May 24 , 1 982 

"Theoretical Analyses of Baroclinic Flows" 
Reporting Period 

August 1 , 1981 - May 31. 1982 

Quarterly Reports Nos. 5-7 
Contract NAS8-33386 . 

Dr. Basil Antar, Principal Investigator 

Prepared for George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


During the present reporting period work was finished on extending 
the Hadley Cell model to investigate the symmetric baroclinic instability 
of the model.’ This work was written in the form of a paper which will 
appear in the June 1982 issue of The Journal of The Atmospheric Sciences. 

A copy of the preprint of this paper is enclosed as Appendix. At the 
present work is being continued to investigate the baroclinic instability 

of an inclined wave for the same Hadley pell model. 

^ 1 ^ . •' * 

.Also during the reporting period work has begun on the coding of a 
program to solve the nonlinear baroclinic instability of an Eady model. 

The basic expression technique have been established and all the theoretical 
background for this work has been completed. It is anticipated the trail 
runs of the program should begin in the coming months. .Initially we will 
investigate the nonlinear evolution of a wave of a single zonal wave number 
and the influence of its hdrmonics. Thenonlinear interaction of more 
than one zonal wave will be investigated subsequently. 
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• Abstract 


A stability analysis of a thin horizontal rotating fluid 
layer which is subjected to arbitrary horizontal and vertical 
temperature gradients is presented. The basic state is a non- 
linear Hadley cell which contains both Ekman and thermal boundary 
layers; it is given in closed form. The stability analysis is 
based on the linearized Navier-Stokes equations, and zonally 
symmetric perturbations in the form of waves propagating in the 
meridional direction are considered. Numerical methods were used 
for the stability .problem. The objective of this investigation 
was to extend previous work on symmetric baroclinic instability 
with a more realistic model. Hence, the study deals with flows 
.for which the Richardson number (based on temperature, and flow 
gradients at mid-depth) is of order unity and less. The com- 
putations cover ranges of Prandtl number 0.2 ^ a <_ 5, Rossby 

“2 2 —4 -1 

number 10 ^ e ^10 and Ekman number 10 -^E <10 .It was found, in 

agreement with previous work, that the instability sets in when 
■the Richardson number is close to unity, and that the critical 
Richardson number is a non-monotonic function of the Prandtl 
number. Further, it was found’ that the critical Richardson number 
decreases with increasing Ekman number, until a critical value 
of the Ekman number is reached beyond which the fluid is stable. 

The principal of exchange of stability was not assumed and growth, 
rates were calculated. A wavelength of maximum growth rate was 
found. For our model overstability was not found. Some com'- 
putations were performed for Richardson numbers less than zero. 


No discontinuities in growth rates are noticeable when the 
Richardson number changes sign. This result indicates a 
smooth transition from symmetric baroclinic instability to a 
convective instability. 






1, Introduction 

It is well known that statically stable baroclinic flow 
can be destablished by two different mechanisms. The first 
occurs in what is usually known as ordinary baroclinic in- 
stability in which the perturbations of maximum growth rate 
have a wavelength in the zonal direction and the gravest 
possible structure in the meridional direction. (Charney, 

1947; Eady, 1949) . This instability is well understood and 
has been observed experimentally (Hide and Mason, 1975) . The 
second occurs in what is known as symmetric baroclinic in- 
stability. In this case, the instability mechanism is such 
that the perturbations of maximum growth rate have meridional 
structure but no zonal structure. This study is concerned with' 
the latter form of instability. 

Two of the earliest analyses of symmetric baroclinic instability 
are those of Solberg (1936) and Kuo (1956) . However, the first 
definitive work is that of Stone. In a sequence of papers, 

Stone (1966, 1970, 1971) examined the stability of a parallel 
baroclinic flow with respect to infinitesimal, three-dimensional 
perturbations . Constant vertical shear and constant temperature 
gradients were assumed and thermal and viscous diffusion effects 
were neglected. Stone found that when the Richardson number, 

Ri, is between 0.95 and 0.25, waves whose axes are in the zonal 
direction have the maximum growth rates. Stone et al . (1969) 
and Hadlock et al. (1972) attempted to observe symmetric baroclinic 
instability, in the laboratory, and some evidence of the predicted 


lneridionf^X structure was seen but a comprehensive experimental 
study of this instability has yet to be performed. 

McIntyre (1970) performed an analysis of symmetric 
baroclinic instability which included both vertical and 
horizontal shears and allowed for the Prandtl number/ fl. The 
model considered was infinite in both the horizontal and vertical 
directions. McIntyre found, in the^ absence of horizontal shear, 
that the maximum value of Ri at which symmetric instability sets 
in (the critical Richardson number, Ri^) is a function of c and 
can be larger than 0.95. McIntyre also found that the most un- 
stable wavelength is infinite and t)iat both monotonic and 
oscillatory instability can occur. In a related analysis, 

Walton (1975) considered only vertical shear in a finite depth 
fluid. When weak viscous effects are included, Walton found 
that Ri_ is lower for both the monotonic and oscillatory ’ in- 

w 

stabilities. He also found that the most rapidly growing modes 
display a weak functional dependence on the diffusion coefficients 
when those are asymptotically small. 

Recently, Emanuel (197’9) presented a detailed study of 
the li ^ar stability of symmetric perturbations in a baroclinic 
flow with both horizontal and vertical shears. Retaining both 
the viscous and thermal diffusivities for a vertically bounded 
flow, Emanuel studied: i) hydrostatic disturbances and ii) non- 

hydrostatic disturbances in a neutrally stratified fluid. He 
* 

found that for both cases the critical Richardson number depended 
on both the thermal and the viscous diffusivities. He also 
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# 


established the fact that the most unstable normal mode is 
determined primarily by the depth of the unstable domain and 
the slope of the isentropic surfaces rather than by the 
diffusive properties of the fluid. Weber (1980) in a some- 
what similar ^^nalysis, examined the symmetric instability of 
a basic state with a constant vertical shear in the presence 
of horizontal boundaries while taking full account of the 
viscous effects for arbitrary ' vertical (stable) stratification. 
Although, he made the simplifying assumption of replacing the 
vertical viscous and diffusive tarms in the perturbation 
equations by a wave mode, his results compare favorably with 
those of Emanuel. It should be noted that in both of the above 

studies the analyses were restrictive; Emanuel investigated only 

# * * 

neutral modes by invoking the principal of exchange of stability 
and Weber used an approximate solution to the perturbation 
equations. Emanuel's paper contains an excellent survey of 
previous work on this problem. 

Previous workers have postulated' syitunetric baroclinic 
instability as the cause of certain physical phenomena. Stone 
(1967) conjectured that the banded structure of Jupiter's 
atmosphere is due to, this instability. Bennets and Hoskins 
(1979) attributed the origin of rain bands and squall lines to 
this instability. Weber (1980) suggested that the generation of 
roll vortices in the oceans and the atmosphere is due to symmetric 
baroclinic instability. 
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i. 

In all of the theoretical studies cited above, the basic 
state whosij stability was investigated was a postulated state 
in which both the shears and the temperature gradients in the 
vertical and the horizontal directions were taken to be constants 
throughout the fluid. In an experimental apparatus whose upper 
and lower boundaries are stationary, such a basic state is a 
good approximation in the interior of the fluid only; it is not 

>r, 

a good approximation near the boundaries, especially when the 
Ekman number is not vanishingly small. In the study, presented, 
in this paper, a fluid contained between two horizontal plates 
of infinite extent is considered and the influence of the Ekman 
and the thermal layers on symmetric baroclinic instability is 
investigated. The analysis uses a realistic basic state which is 
obtained through the solution of the governing nonlinear equations of 
motion for abitrary external forcings, e.g., the Ekman number, 
the vertical stratification, etc. This basic state was used 
in a recent paper by the authors in an investigation of ordinary 
baroclinic instability (see Antar and Fowlis, 1981) . 

It is obvious that if a direct analog of a laboratory ex- 
periment is required, then this model .is not sufficient since it 
does not allow for the side walls and the side layers. In a 
recent 'study Quon (1980, 1981) simulated laboratory conditions 
for the onset of symmetric baroclinic instability through the 
numerical solution of the full Navier-Stokes , energy equations 
for .the flow in a cylindrical differentially heated, rotating 
annulus. Although such an analysis is of great value in. modeling 
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the true experimental conditions, the amount of computer time 
required for the calculation of each case precludes a thorough 
probing of the parameter space for the criteria for the onset 
of the instability. It is our opinion that the present study 
fills a needed gap between the constant shear, constant tem- 
perature gradients analyses discussed earlier and that of the 
exact solution of Quon. 

The basic state model is reviewed in Section 2 . The 
perturbation equations for the symmetric instability and their 
method of solution are presented in Section 3. In Section 4 
we discuss the results of the present analysis and 
how they compare with the previous work. Our conclusions 
are presented in Section 5. 




2 , The Basic State 

The basic state used in the present stability analysis 
is identical to that used in the work of Ax^itar and Fowlis (1.981) 
(hereinafter referred to as I) . This model is described only 
briefly here; for further details and justifications, the 
reader is referred to I. We consider a Boussinesq fluid confined 
between two horizontal plates which are set a distance, d, apart. 
The coordinate system used is rectangular Cartesian with axes 
(x,y,z) corresponding to the eastward, northward and vertical ‘ 
directions, respectively, and with the oi’igin midway between the 
plates. A sketch of the model is shown in Fig. 1. The plates 
are assumed to extend to infinity in the x-directipn and to 
large distances in the y-direction. The plates and the fluid 
•are taken to rotate as a whole about the vertical axis with a 
constant angular velocity, n. A tfl?mperature gradient in which 
the temperature decreases linearly in the y-direction is main- 
tained on both plates. To assure that the vertical stratification 
and hence the Richardson number, is arbitrary and externally 
fixed, the temperatures of the- upper and lower piates are set to 
. differ uniformly by a constant amount AT for all (x,y). 

The velocity and temperature fields are governed by the 
Navier-Stokes , energy and mass conservation equations. Under the 
conditions assumed here for the basic, two-dimensional, steady 
state, these equations, in a rotating reference frame, reduce to 
the following set: 

• 6 


( 1 ) 




-2V « E3^U/B25=^ , 

2U * E3^V/3z^ - 3p/8yr (2) 

V3T/3y = (E/ae) 3*^T/3z^/ (3) 

T = 3p/3z, (4) 

where is the velocity vector (U/V) , p is the pressure, T is 
the temperature and 
e = 

E *= , 

a « v/k, 

are a thermal Rossby number, the Ekman number and the Prandtl 

number, respectively. In Eqs. (l)-(4), length, time, velocity 

' ■ . ■ -1 • 

and. temperaturfe have been made dimensionless usrng d, n , 

a gd y/n, and where a is the coefficient of thermal expansion, 
y is the imposed horizontal temperature gradient and g is gravity. 
The main assumption's used in arriving at Eqs. (l)-(4) are that 
the vertical velocity component, W, is negligible and that 
3T/3y = -y throughout the region of interest, (see I). 

The solution of the system (l)-(4) which is consistent 
with the no s3>ip and perfectly conducting conditions at the solid 
boundaries, i.e. , 

U*V=0, T=+ AT/2-y 0 z = + 1/2, (5) 

is given by 

•U (z) = - f (z)/8 + z/2 , ■ (6) 
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V(z) « g{z)/8, (7) 

T(y,z) » -y zAT + ea[2z - f(z)J/8, (8) 

where 

f(z) « [coshR(z+l/2) cosR(z-l/2) - coshR(z-l/2)cosR(z+l/2) ] /h (R) / 

g(z) = IsinhR(z4-l/2) sinR(z-,l/2) - sinhR(z-l/2)sinR(z+l/2) ]/h (R) / 

h(R) »= sihh*(R/2) + sin^(R/2)/ 

R = 

The velocity and temperature profiles given by (6) -(8) are 
shown in Pig. 2 for different values E, c, a, and AT. The 
limits of applicability of this model are discussed In I. 
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3. The Perturbation Equations 

The instability problem is tackled by imposing perturba- 
tions on the basic state presented above. Since we are concerned 
with symmetric perturbations, the perturbation functions will be 
proportional to exp [ii(y ” ct) ] which represents a wave with 
meridional wavenumber propagating in the meridional direction with 
speed c. Introducing the perturbations into the Navier-Stokes , 
energy ahd mass conservation 'equations and upon linearizing the 
perturbation equations we obtain the following equations for the 


perturbation functions , 

E(D^-Jl^)^w + ijl[(c-eV) (D^“£^)w + ewD*V] -2iADu-A^Q = 0, (9) 

E(D=^-jl*) u + iJl (c-cV) u-ewDU - ~ Dw = 0, (10) 

E(D*-Jl*) e + iS.a (c-eV) e- -—cDw - ecwDT =0, (11) 

i£u‘ + Dw = 0, (12) 

where 
D = d/dz. 


In the above the perturbation velocity is = (u,v‘,w) and 6 
is -the perturbation temperature. Note that the perturbation 
equations are not hydrostatic . 

To complete the specification of the problem# the following 
boundary conditions a.re imposed on the perturbation functions. 

u=w=Dw=e=0@z=+ 1/2. (13) 
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The perturbation problem as posed constitutes an eigen- 
value problem. We take 1 to be real and c * c + ic. 

£ 

to be complex^ where c denotes the propagation speed and c., 

* X 

the growth rate. Monotonlc instability occurs for ^ 0 while 
Cj. “ 0, whereas overstability occurs for c^ ■ 0 while 0. 

It is tempting at this stage to invoke the principle of exchange 
of stability and to set * 0 hence look for neutral 

modes alone and reducS the number of parameters involved. Emanuel 
(1979) adopted this procedure. However, since McIntyre (1970) 
had established that for his model (see Section 1) oscillatory 
instability exists, it was decided to leave and as eigen- 
values of the problem. Using this method the analysis will allow 
us to compute growth rater and to determine if overstability 
exists . . . . ■ 

The solution to the problem (9) - (13) is obtained through 
numerical integration of the full eighth-order differential 
system. This procedure was necessary because the coefficients 
in the equations (i.e., the basic state) are functions of the* 
height, If we wish to account fully for the Ukman and thermal 
layeifs, the solution cannot be obtained in a closed form. Note, 
that if the basic state profiles had been taken to be linear in z, 
as was done in most of the previous work on the symmetric baroclinic 
instability problem, then the set of governing equations would 
reduce to a differential system with constant coefficients whose 
solution could. be determined in a closed form. 
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The numerical solution was obtained using the shooting 
technique. The code was identical to the one used in and 
more details of the numerical work are given in I. Briefly, the 
shooting technique was enhanced with the orthonormalization 

I 

method, allowing the stability analysis to .be performed even 
when the Ekman number is very small. A Newton -Raph son method 
was implemented for convergence on the eigenvalues, and an eighth- 
order Runge-Kutta method was used for the integration of the 
initial value problem. 


4. 


Results and Discussion 


The results which follow constitute the solution of the 
eigenvalue problem posed by equations (9) - (11) with the 
boundary conditions (13) for prescribed values of the 
parameters. The parameters of the problem which appear ex- 
plicitly in the equation's are, e, the thermal Rossby number 
E, the EJcman number, a, the Prandtl number, AT, the imposed 
vertical temperature difference, and £ and c, the wavenumber 
and wave speed, respectively. Note that the first four 
parameters appear in the basic state as well as the perturbation 
equations, but the last two belong to the stability problem 
alone . 

As was discussed in Section 1, the Richardson number plays 
a significant role in the symmetric baroclinic instability problem. 
While this number does not appear explicitly in our formulation 
of the basic state, it may be defined using the parameters which 
do appear explicitly. The Richardson number is defined by 


Ri = 


gg 3T*/3z* 


(14) 


where 3T*/3z* and 3U*/3z* are the* dimensional basic state tem- 
perature and zonal flow gradients, respectively. It is obvious 
from this definition and from expressions (6) - (8) that for the 
present model Ri is not a constant but is a function of z. We * 
may define a "bulk" Richardson number by the value of Ri at 
mid-depth between the horizontal plates, viz. 


P,* ~ rvrr 9T*/3 z* 


Z K 0 • 


In all of the subsequent discussion, Rij^ will be used instead 
of Ri and the subscript b will be dropped. We shall see that 
this definition is not only convenient, but it has dynamical sig- 
nificance in this problem. Substituting the temperature and 
velocity gradients for the basic state, expressions (6) and 
(8) , respectively, into the above definition we obtain^ 



,AT + ae. 
' 4 ^ 


(15) 


which is a constant and defined exclusively by the basic state. 

The above definition shows that Ri is a function of the 
explicit parameters of .* the problem.' This relationship was’ used 

* » • f 

when choosing values of e, AT and a for each eigenvalue search 
in order to be sure that -we were within the correct range of Ri 
for instability. Stone (1970) and McIntyre (1970) showed that 
symmetric baroclinic instability occurs for Ri close to unity 
(see Section 1) . The above definition indicates that for a 
given value of c, there is a range of. e and of AT for which -the 
critical value of Ri is achieved. 


For (j-1, expression (15) specifies that AT must be less 
than zero to obtain the desired range of values of Ri for in- 
stability regardless of the value of e. However, setting AT< 0 
arbitrarily could lead to an unstable vertical stratification 
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in the interior and hence a negative value of Ri. Instabilities 
occuring under the above conditions with negative Ri are not 
classed as symmetric baroclinic instabilities and are outside 
the range of interest of this paper. Such instabilities 
have been discussed by Hathaway et al, (1979/1980), To obtain symmet 
baroclinic instability for 0*1/ AT must be set negative with the 
additional constraint of |4 AT/e| <^1». Clearly/ there exists 
a range of values AT which satisfies these two criteria. We decided 
to proceed with AT « -0.2 for o * 1. Fig. 3 shows that this 
value of AT keeps the interior vertical stratification stable. 

The above arguments reveal that the correct choice of AT for 
symmetric instability is a delicate matter and depends strongly 
on 0 . Expression (15) should be of great value in guiding ex- 
perimental searches for symmetric baroclinic instability in the 
laboratory. 

Pig. 4 shows the growth rates for different values of the 

meridional wavenumber/ li as a function of e for E * 10 / AT = 

-0.2 and a * 1. Note the existence of unstable modes over a 

finite range of e. This feature resembles the growth rate 

properties of ordinary baroclinic instability. The unstable 

region in e is bounded by two stable regions/ one in which Ri is 
• » 

close to unity and the other in which Ri is small or negative 
(Ri * 1 - O.B/e) . 

Fig. 4 also shows that as c varies, the growth rate has a 
maximum value for a specific value of the wavenumber. Fig. 5, 
which is for the same values of the parameters as the results in 


Pig. A, is a plot of the maximum growth rate for each wavenumber 
versus wavenumber. The greatest growth rate occurs for about 
i * 10. This result is at variance with the findings of Stone 
(1970) , who found that the maximum growth rate is achieved as 
i w, and McIntyre (1970) , who found that it is achieved as 
i ■+• 0. However, since both of these previous studies are for 
essentially inviscid flow, it is concluded that including the 
viscous effects determines the wavelength of maximum growth rate 
for the symmetric baroclinic instability problem. 

Figure 6, which is for the same values of the parameters as 
the results in Pigs. 4 and 5, shows neutral stability curves in 
the (E - e) number plane for three different wavenumbers. A sig- 
nificant feature of these curves is the existence of the maximum 
(critical) value of E‘ above which ail modes are stable. Another 
feature is the existence of two stable regions, or two branches of 
the neutral stability curves, for all values of E smaller than the 
mode critical values. Note that these curves are loci of points 
for which c^ *= 0. For all the modes shown in this figure -and 

I 

those in the previous and subsequent figures, the value of the wave 
speed, c^, was always found to be zero, indicating that these are 
monotonic instability modes and that the principle of exchange of 
stability holds for this problem. .Throughout this work we did 
not find any oscillatory unstable modes. Thus, it is doubtful if 
overstability exists for the present problem. 

The curves in the upper portion of Pig. 6 indicate that 
e “ as E -► 0 for all of the wavenumbers shown. To observe this 
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criterion in terms of Ri, the neutral curves of Fig. 6 are plotted 
again in the (E-Ri) plane in Fig. 7. Pig. 7 shows that as E 0, 
the upper branches of the curves approach Ri «* 1 for all wave- 
lengths and that the region of instability increases with in- 
creasing wavenumber. Weber (1980) / in his model for symmetric 

I 

baroclinic instability (see Section 1 ) , found the critical value 
of Ri, Ri , to be a monotonically decreasing function of E as the 

w 

latter increases from zero. The results in Fig. 7 show that this 

fundamental dependence is correct up to the value of E^, but above 

that value there is a region of absolute stability regardless of 

the value of Ri. We believe that the existence of E is due to 

c 

the accurate inclusion of the viscous* and thermal dxffusive effects 
in our model. None of the previous work included both of 
these effects properly. 

To investigate the effects of the Prandtl number on the 
stability criteria discussed above, we proceeded to find solutions 
for 0 *= 5. As was discussed earlier, changing the value of a 
necessitates a search for the appropriate value of the vertical 
temperature difference, AT, such that Ri is in the proper range 
for symmetric baroclinic instability,. . The value chosen for o *= 5 
was AT = -4.5, which allows Ri to be in the vicinity of unity 
(Ri = 5 - 18/G). Fig. 8 'shows the growth rates as a function of 

-3 

e for the above values of 0 and AT and for E = 10 . Again, it 

is seen that the unstable modes separate two stables regions 
in e (or Ri) . Fig. 8 shows that the maximum growth rate increases 
with increasing wavenumber, A. In fact, the total range, of 
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instability increases with Note that the upper bound of 
instability in e first increases with increasing wave numbers 
up to A. * 4 and then decreases, reaching an asymptotic value of 
e » 4.55 for X. > 6 which corresponds to Ri * 1 for the parameters 
used. Another feature# which is different from the o « 1 case, 
is that the maximum c*' 'th rates occur for approximately the same 
value of e for f .. cne wavenumbers shown. For this case this 
value of e translates to a value of Ri » 41 

Pig. 9, which is for the same values of the parameters as 

the results in Fig. 8, shov;s the variation of the maximum growth 

rates as a function of the wavenumber, H. This figure, reveals 
that there is still a wavenumber of greatest growth rate which 
occurs about % » 13.5. However, the decrease in the growth 
rate is much smaller for wavenumbers above the maximum value 

than it is for the a * 1 results (see Fig. 5) . 

The neutral stability curves in the (E - Ri) plr.nc are 

shown in Fig. 10 for X. = 1 and 4. In this figure only the positive 

part of the Ri range is shown. Note that there are two values of 

Ri as E 0, each corresponding to one wavenumber, and that the 
c 

range of instability is much larger for ^. = 4 than '‘■'hat for X- * 1 . 

To appreciate more fully the influence of the Prandtl 
number on the symmetric instabiity criteria. Pig. 11 which shows 
the neutral stability curves for X = 4 in the (E - Ri) plane for 
a » 1 ( AT « -0.2) ,0=2 (AT = -1/5) and o - 5 (AT = -4.5) was. 
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prepared. Except for a = 1, only the upper branches of the 

curves appear for Ri > 0. All the curves shown tend to specific 

asyiaptotic values of Ri as E 0. For o “ 1, the asymptotic 

o 

value of Ri_ is close to 1# while for a *= 2 and 5, the asymptotic 

values are close to 1.075 and 1.695 respectively. The latter 

two values are significantly greater than !•. Thus, the trend 

for the asymptotic value of Ri is co increase with increasing 

0 for 0 > 1. This fundamental depeiu\tnc{.- is in agreement with 

work of McIntyre. fl970), although his model is somewhat different 

from the model being considered here (See Section 1) . McIntyre 

found that the asymptotic value of Ri is given by 

c 



(1 + 0 ) 
40 


2 


(16) 


• For 0*1, 2, and 5, this relationship yields values of Ri . of 1,' 

... * * w * ' 

1.125 and 1 . 8, respectively . These values are close to the values 
found from Fig. 11. 

Fig. 11 also shows that a;s 0 increases, the whole region 

of instability is markedly increased. The value of E increases 

c 

drastically as 0 is increased but there is no observable difference 

in the values of E for a = 2 and 5. . . 

c 

Next, we present results for 0 < 1. Again, the values for 
AT were chosen such that the range of Ri includes unity and values 
less than unity. Expression (15) indicates that when 0 < 1 and 
AT <■ 0, the maximum possible value for Ri is 0 and thu«ij .;he range 
does riot include unity. However, when o < 1 and AT >0, values of 
Ri can be obtained which are both greater and less than unity,' We 
proceeded with AT > 0 . 
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Fig. 12 shows the neutral stability curves in the (E - Ri) 

plane for a = 0,5 (AT »= 0.2) and a * 0.2 (AT >= 0.5). The curves 

are for A “ 12, the value which yielded the maximum value for 

Ri i^hen all the other parameters were held fixed. Again, the 

values for Ri are seen to increase to specific asymptotic values 

as E -► 0. Note that the tendency to an asymptote is not as evident 

for a *= 0.5 as it is for a * 0.2.-' T?he main features of Fig. 12 

are the decrease in the value of Ri and the decrease in the ‘total 

c 

instability area as a increases to unity. The results also in- 
dicate an increase in E^ with increasing cr. 

-4 

Fig. 12 gives the values of Ri^ at E = 10 for a * 0.5 

and 0.2 as 0.946 and 1.21, respectively. Pig. 13 is a comparison 

-4 

of all our results for the largest asymptotic (E * 10 )• values 

of Ri with McIntyre's results as a function of a. The dots are 
o 

our results and the continuous curve is a plot of expression (16) . 

Note that McIntyre's criterion yields a symmetric dependence of 

Ri as a function of a about a = 1. The results of the present 
c 

study agree well with McIntyre for cr ^ 1 but give lower values 

for Ri for a < 1. 

c • • 
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5. Conclusions 

The stability of a nonlinear rotating hci .ley cell has 

been examined with respect to symmetric perturbations. This 

Hadley cell was derived previously (see I) and contains both 

Ekman and thermal boundary layers. The stability analysis is 

based on the Navier-Stokes equations and considers infinitesimal 

waves with vertical and meridional structure but no zonal structure. 

The stability analysis was performed numerically. For examining 

the parametrip behavior of symmetric baroclinic instability, the 

work is an improvement over all previous work since a realistic 

basic state is used and since the stability analysis includes 

'both viscous and thermal diffusion effects. 

The computations cover the following ranges of the non- 

“2 2 —4 “I 

dimensional parameters*-. 0.2 ^ a <5, 10 e <..10 / ,10 <. ^ 1 10 f‘ 

-4 . 5 ^ AT <_ 0.5, and -100 £ Ri £ 2. It was found that the instability 

sets in when is close to unity and that the critical Richardson 

number, Ri , is a strong function of both a .and E for the cases 
C * 

considered. We extra^fxslated from these cases to obtain values of 
Ri for the inviscid limit, E *► 0 , and found that Ri_ has a maximum 
value close to unity for o = 1,’ -The general trend of these results 
is in agreeme.’-it with McIntyre (1970) and for £ 1 the agreement is 
very good but fora < 1 it is not so ’good. We should not be suprised 

to find differences between the present and previous work because 
of the differences in the basic states and other aspects of the • 
models'. 


We also found for fixed a that Ri„ decreases with in- 

c 

creasing E until a critical value of E is reached beyond which 
the fluid is stable. The principle of exchange of stability was 
not assumed and growth rates were calculated. A wavelength of 
maximum growth rate occunring for small wavelengths was found. 

For our model overstability was not found. Some computations were 
performed for Ri < 0. In particular, for a = 5 the growth rate 
maximum occurs at about Ri = -4 . No discontinuities in the growth 
rates are noticeable when Ri changes sign. This result indicates 
a smooth transition, from symmetric baroclinic instability to a con- 
vective instability. 
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Legend 

1. A sketch of the geometry of the model and the 
coordihate system. 

2a. The zonal velocity component of the basic state ^ 

U/ as a function of z, for three values of E. 

2b. The meridional component of the basic state, V, 
as a function of z, for the same values of E as 
the results shown in Figure 2a. 

2c. Basic state temperature profiles, T at y * 0 for 
the same values of E as the results shown in 
Fig . 2a and for c = 1, e, = 1 and AT «= 0. 

3. Basic state temperature profiles, T, at y = 0 for 

three values of AT and for c = 1, e = 1, and 

E »= 5 K . 10 ^ , • 

4. Growth rates, Ic^, as a function of e for different 

-3 

values of -i and for a = 1, AT * -0.2 and E = 10 

The upper abscissa shows the values of Ri corresponding 

to E . 

5. The maximum growth rates as a function of the 
meridional wave number, i, for the same conditions 
as Fig. 4. 

6. Neutral stability curves, c^ = 0, as functions of 

E and e for different values of SL and for a = 1 and 
AT = -0.2. 

7. The neutral stability curves of Fig. 6 as functions 
of E and Ri. 
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Figure 8. 


Figure 9. 

% 

Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Growth rates, as a function of e for different 

values of i and for o « 5, AT * -4.5 and E * 10“^. 

The upper abscissa shows the values of Ri corresponding 
to e. 

The maximuin growth rates as a function of A for the 
same conditions as Fig. 8. 

Neutral stability curves as functions of E and Ri 
for different values of t and for a ® 5 and AT » -4.5. 
Neutral stability curves as functions of E and Ri for 
different values of a (^1) and for i = 4. The 
dashed lines are McIntyre's results. 

Neutral stability curves as functions of E and Ri for 

» 

different values of a (< 1) and for a = 12. The 

dashed line is McIntyre's result, 

A comparison between the critical value of the 

Richardson number, Ri_/ as calculated from Eq. (16) 

c 

( -) and the results presented in this paper 

-4 

as a function o.f a for E =. 10 (O). 
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